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SUMMARY  AND  PREFACE 


The  X  and  Y  functions  of  Chandrasekhar  are  of  great  importance 
in  the  theory  of  multiple  scattering  in  a  finite  slab.  Their  properties 
are  best  determined  analytically  through  the  use  of  integral  equations. 
From  the  computational  view,  though,  they  are  best  treated  as  solutions 
of  a  system  of  integro-dif f erential  equations. 

Our  aim  is  to  produce  a  derivation  from  first  principles  which 
illuminates  the  nature  of  the  integro-dif ferential  equations  and  which 
will  aid  us  in  our  computational  studies  of  transport  processes.  A 
knowledge  of  ir  2gral  equations  and  Chandrasekhar's  book,  Radiative 
Transfer , ^  is  assumed. 
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I.  INTRODUCTION 

(2) 

In  a  preceding  paper, '  '  with  the  aid  of  the  invariant  imbedding 
technique,  we  derived  the  functional  equations  for  the  Fredholm  resol¬ 
vent  for  the  photor.  diffusion  equation  in  the  case  of  the  slab  geometry. 
In  the  present  paper,  making  use  of  the  relationship  between  the  source 
function  and  the  resolvent  obtained  previously,  we  show  a  new  approach 
for  finding  the  integro-diff  ere..cial  equations  for  the  X  and  Y  functions 
of  Chandrasekhar^  in  the  theory  of  radiative  transfer  in  a  finite 
flat  layer.  Since  these  equations  are  useful  for  the  numerical  cal¬ 
culation  of  the  X  and  Y  functions,  they  play  a  significant  role  in  the 
inverse  problem  of  light  scattering  in  slab  geometry. ^  Let  us  also 
note  that  there  is  no  problem  of  uniqueness  in  the  case  of  conservative 
scattering  when  vie  employ  the  integro-differential  equations  for  the  X 
and  Y  functions  rather  than  the  integral  equations. 

II.  MILNE*  S  INTEGRAL  EQUATION 

Consider  the  diffuse  reflection  and  transmission  of  light  by  a 
plane-  arallel,  homogeneous,  and  isotropically  scattering  atmosphere 
of  the  finite  optical  thickness  Tj  >  0  with  the  distribution  B^(t) 
of  internal  emission  sources.  Let  a  parallel  beam  of  radiation  of 
constant  net  flux  ttF  per  unit  area  normal  to  the  direction  of 
propagation  be  incident  on  the  surface  t  *  0  in  a  fixed  angle 
cos"Vq(0  <  iiQ  s  1)  with  the  inwards  normal. 

Following  the  notation  of  Chandrasekhar,^  we  shall  use 
I(t,+  h)  (0  <  M  1)  t0  denote  the  intensity  of  radiation  at  level 
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In  the  direction  +  p  directed  towards  the  top  surface  t  *  0,  and  further 
I(T>-|i)(0  <[j,  s  1)  for  that  directed  towards  the  bottom  surface  t  - 

The  equation  of  transfer  in  the  diffuse  radiation  field  takes  the 

form 

+1 

p,  \  \  J  I(T»M>#)d^ ' "  B(t)  ,  (1/ 

where  \  is  the  albedo  for  single  scattering  and  B(t)  is 

B(t)  -  BjCt)  +  |x  F  e'T/^  •  <2) 

Equation  (1)  should  be  solved  subject  to  the  boundary  conditions 

I(Of-  p.)  “  0  and  1(t^,+  M>)  *  0  (0  <  (J>  £  1)  •  (3) 

Equation  (3)  represents  the  fact  that  no  diffuse  radiation  is 

incident  on  the  surfaces  from  the  outside. 

Starting  from  the  formal  solution  of  Eq.  (1)»  we  have  the  Milne 

integral  equation  for  the  source  function  S(t) 

s(t)  -  T :  \{S(t)J  +  B(t)  ,  <4> 

where  the  truncated  Hopf  operator  is  given  by 

T1 

*T{f(t))  -  |J  f(t)E1(|t-Tj)dt  (0  <  T  <  TX)  • 

T  0 

In  Eq.  (5)  Ex  is  the  first  exponential  integral 

Kj(T)  -  J  e-T/“  *  . 

0 


(5) 


Alternatively,  Eq.  (4)  is  written  in  the  form 


X 

S(t)-B(t)+J  K(T,t;T  ,  (6 

0 

where  the  Fredholm  resolvent  K(T,t;Tl)dTdt  is  the  probability  that  a 
photon  emitted  in  the  optical  depth  interval  (t,t+dt)  will  be  re¬ 
emitted  in  the  optical  depth  interval  (T,T+dT)  after  one  or  more 
scattering  processes. 

Then,  once  the  resolvent  K  has  been  determined,  Eq.  (6)  enables 
us  to  evaluate  the  source  function  which  gives  us  the  intensities  of 
the  emergent  diffuse  radiation  as  follows: 

T1 

K0.+ 10  -  J  3(0*-*'“  £  ,  (7) 


I(t1#-  M>)  -  J  S(t)e —  . 


The  directly  transmitted  intensity  is  |  e  6(u-u  ) 

2  * 


FREDHOIM  RESOLVENT  OF  THE  PHOTON  DIFFUSION  EQUATION 


/  j} 

In  a  preceding  paper,  with  the  aid  of  the  invariant  imbedding 
technique,  we  obtained  the  following  equations  (making  use  of  the 
optical  homogeneous  properties  of  the  medium): 
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Then, 


K(x,y;Tl)  -  K(j:-y,0;r1)  +  J  ^K(t+x-y,0;Tl)K(0,t; 

0 

-  K(T1-t-x40S0;Tl)K(0,T1-t;T1)jdt  , 


where  x  >  y. 

With  the  aid  of  Eq.  (11),  the  resolvent  K(x,y;r^)  is  expressed 
in  terms  of  K(x,0;t^),  a  significant  reduction.  Equations  (9),  (10), 
and  (11)  reduce  to  those  given  by  Sobolev.  ^ 

On  the  other  hand,  by  means  of  another  analytical  procedure, 
we  know  that 

T1 

K(x,y;Tp  -  k(x,y)  +  J  K(x,y /;x1)k(y ',y)dy '  (12: 


where 


k(y  \y)  »  2  ^  Vly'-y|>  * 


Furthermore,  from  Eq.  (4)  with  B^t)  *  0  and  F  -  4,  we  have1 

v 

*(x,T1)  -  g(x)  +  J  *(x#,T1)k(x/,x)dx/  , 


where  k(x,x')  is  given  by  Eq.  (13),  and 


«(x,T1)  -  “  J  S*(x,&l)  ^  , 


gOO  -  \  \  Ex(x) 


(16) 
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In  Eq.  (15)  the  source  function  S  (x,n)  satisfies  the  auxiliary 
equation 

a?) 

where  1  on  the  left*hand  side  represents  an  identity  operator.  The 
* 

quantity  S  (t,p.)  corresponds  to  the  probability  that  a  photon  absorbed 
at  level  t  will  reappear  in  the  direction  +  ^(0  <  y,  s  1)  in  the  ra¬ 
diation  emerging  from  the  surface  t  =  0.  Equation  (17)  is  also  called 

a  photon  diffusion  equation. 

(2) 

In  a  preceding  paper,  it  is  also  shown  that 

K(x,0;Ti)  *  $(x,Tl)  ,  K(t1-x,0;t1)  -  SCj-x.Tp  .  (18) 
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Putting 

B^t)  -  0  and  F  -  4  ,  (21) 

Eq.  (20)  for  S  (T,pi)  becomes 


as  (t,h) 

~"^7~ 


X 

2  S  (Tl>^  J  S  (T1”T,n')  , 


(22) 


which  reduces  to  that  given  by  Sobolev,  ^  Busbridge'' >J  and  Ueno 
Writing 

S*(0,^)  -  i  x(n,Tl),  S*(T : ,&i)  =  I  YCi,Tl)  , 


.(7) 


(8) 


(23) 


and  letting  t  *  0,  from  Eq.  (22)  we  find  the  desired  equation 

) 

Equation  (24)  is  one  of  the  requisite  integro-differnntial 
equations  for  the  X  and  Y  functions. 

Now  we  derive  the  other-  Let  the  Milne  integral  equation  for 
the  source  function  S(TpT)  be 

S(T]“t)  -  a  A^{S(T;l- t)  +  B(t^-t)  .  (25) 

With  the  aid  of  the  Fredholm  resolvent,  Eq.  (25)  is  written 
in  the  form 

T1 

S^-t)  =  B^-t)  +  J  K(T,t;T1)B(T1-t)dt.  (26) 

0 

On  differentiating  Eq.  (26)  with  respect  to  we  get 


axOi,T  )  1 

-jT— -5T:yoi.Ti)  Jtoi'.Tj 

o 
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dS(r1-T) 

3t, 


^B(t.-t)  p1  3K(T,t;T-) 

-TT—  +  K^VTi>B<°>  +  j  — 5T - 

1  0  1 

T 

dB(T.-t) 


B(T^-t)dt 


(27) 


+  j  at  . 


Making  use  of  Eq.  (21),  from  Eq.  (25)  we  have  the  auxiliary 

★ 

equation  for  the  source  function  S  (t.-t,^) 


*  _  -(t.-t)/h 

[1  ’  X  A]t{S  (Tj-t.iA)}  -  X  e 


(28) 


The  combinations  of  Eqs.  (27)  and  (28)  leads  to 

dS*(T  -T,n)  -(t,-t)/h 

-  -  7  X  e  +  K(t,t1;t1)x 


di\ 


M* 


1  aK(T,t;Tl)  -(Tj-O/H  j  1 


Bt, 


<lt  ‘  a  J  K<T,t;Tj)^ 


-(Tj-O/pi 


0 


dt 


l  r  "Ct,-t)/h  _  pl  -(t  -t)/M. 

*  -  ~  [x  e  +  X  J  K^.tjTj)#  dtj  (29) 

0 

p1  dK(T,tJT  )  -(T.-t)/li 

+  I  [k(t1“t»0»t1)  +  J  — bT -  e  dtJ 

n  1 


-  -J  S*(t1-t,h)  +  \  S*(0*) 


S*(t1“T,H/)  ^7  . 


Then,  letting  t  ■  0  in  Eq.  (29)  and  recalling  Eq,  (23),  we  have 


8Y(h,t,)  .  1  _ 

—srr  ■  -  i  Y(iiiV 


L 

f  ^  .  (30) 
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Equation  (30)  is  the  required  integral  en nation  for  the  X  and 
Y  functions.  Equation  (29)  reduces  to  that  given  by  Busbridge'  * 
and  Ueno,^^  respectively. 
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